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ABSTRACT 

Definitive attitude solutions are supposed to be the most accurate possi- 
ble. For the Earth Radiation Budget Satellite (ERBS), this has been ac- 
complished by using gyro rates to transform many nonsimultaneous 
observations to a common time point and then averaging to reduce the 
effects of observation noise. Rate quality is critical to realizing improved 
accuracy with this method. Gyro deterioration, which shows up as large 
observation residuals and discontinuities between contiguous batch solu- 
tions, now discourages using the batch approach for ERBS. To address 
this problem, a simple Kalman filter is tried in place of the batch estima- 
tor. . The filter works well as long as the attitude is completely observable. 
During periods without Sun coverage, however, the extrapolated yaw may 
diverge and then change abruptly when the Sun returns to the sensor field 
of view. Causes of this behavior are discussed, and some solutions are 
tried that address the observability aspect of the problem. 


*This work was supported by the National Aeronautics and Space Administration (NASA) /Goddard 
Space Flight Center (GSFC). Greenbelt, Maryland, Contract NAS 5-31500. 
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1. INTRODUCTION 


Definitive attitude solutions are to be the most accurate possible and are to be continuous 
over the entire mission. Under the best conditions, this can be difficult both technically 
and operationally (Reference 1). In this paper, we consider one of the technical prob- 
lems— that caused by attitude propagation error. 

As a spacecraft gets old, its gyros may become noisy. This has been the case for the 
Earth Radiation Budget Satellite (ERBS). Without accurate angular velocity measure- 
ments, its definitive attitude determination system has been unable to provide continuous 
solutions or fit all observations. 

To remedy this, the batch estimator is replaced by a Kalman filter that accounts for the 
gyro noise. The results obtained are described, and recommendations are made for fur- 
ther work to improve the approach. 


2. ERBS BACKGROUND 


ERBS is an Earth-pointing spacecraft in a 57-degree (deg) inclination, nearly circular orbit 
with a period of 97 minutes. It has two horizon sensors for measuring geodetic pitch and 
roll plus two Sun sensors whose alpha angles measure yaw. In addition, ERBS has a 
three-axis magnetometer. The horizon sensors are accurate to 0.5 deg, due primari y o 
horizon radiance effects. The Sun sensor is accurate to 0.05 deg, and the magnetometer 
is accurate to 3 deg. Because of its lower accuracy, the magnetometer is not used for 
definitive attitude determination. ERBS also has two redundant three-axis gyros for meas- 
uring its angular velocity. Both gyros have now partially failed. The first failure occurred 
on a pitch-axis gyro in August 1986, 22 months after launch. The second gyro failed on a 
pitch axis as well in July 1988. At present, only roll and yaw rates are available. 


While the horizon sensors always see the Earth, the Sun sensors see the Sun for on y pa 
of each orbit. This means that attitude is not completely observable at all times, and good 
angular velocity measurements are needed to provide a complete history of the attitude. 
This incomplete observability also necessitates special care in dividing up each day s data 
for computing attitude. 


3. DEFINITIVE ATTITUDE DETERMINATION 


The ERBS definitive attitude system has three parts that determine the effect of propaga- 
tion error-the estimator, the segmenter, and the smoother (Reference 2). The current 
batch least squares estimator updates its previous estimate of the state (epoch attitude, 
gyro bias, and scale factor errors) using all the observations at once. Sensor observations 
from different times are in effect propagated back to one epoch time, where they are 
averaged to reduce the effect of measurement noise. 

If propagation were perfect, any number of observations could be so transformed, giving 
a solution of unlimited accuracy. In practice, however, propagation does ad unce ainty, 
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which increases with the time between the epoch and the observation times. Thus the 
residuals, or differences between the observed and predicted values, reflect both sensor 
noise, which is stationary, and propagation error, whose uncertainty grows with time. 

If these residuals grow linearly with time from epoch or in proportion to the gyro output 
rates, they can be reduced by solving for rate bias and scale factor corrections along with 
the epoch attitude. Since bias errors grow with time, knowing the attitude at both ends of 
the batch affords the best opportunity for estimating biases. For this reason, the seg- 
menter divides the day s data so that there is Sun sensor coverage at the start and end of 
every batch. This ensures that the rate biases can be accurately determined. 

The smoother also serves to reduce the evidence, if not the effect, of propagation error 
Because not all gyro errors are eliminated by bias and scale factor corrections, the time 
dependence of the attitude history may still be incorrect. No matter what the choice of 
epoch attitude, the solution cannot be correct over the entire batch. This shows up in the 
definitive attitude as discontinuities between batches. The smoother discards points 
around the junction and replaces them with values obtained by linearly interpolating be- 
tween the new end points as in Figure 1. This is done without regard to sensor observa- 
tions, and may not significantly improve the solution. 

The batch estimator, by ignoring this random part of propagation error, creates a problem 
that the segmenter cannot fix and the smoother simply covers up. 

4. MODELING PROPAGATION ERROR 

Although propagation error can cause batch junction discontinuities, it remains to be 
shown that the ERBS gyro deterioration actually produces such discontinuities. To make 
this connection, a simplified Farrenkopf gyro noise model is used to predict discontinuity 
size based on the observations of gyro noise. These predictions are then shown to track 
the daily average pitch, roll, and yaw junction discontinuities. 

Unlike the white noise seen in the sensor observations, errors in the propagated attitude 
tend to change continuously. For this reason, the propagated angle (0) and its error 
(A0) are modeled as a random process. The model in Figure 2 is a linear second-order 
system with four zero mean value sources of error: 

• Float torque derivative noise (n u ) 

• Float torque noise (n v ) 

• Initial attitude error (A0(O)) 

• Drift rate bias (b) 

The float torque derivative noise is integrated to give a random walk that is added to the 
true angular rate (<y). A constant, but imprecisely known bias and the float torque noise 
are also added at that point. At the output, another constant error is added representing 
the initial or epoch attitude error. * 
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Figure 1. Smoothed Junction Discontinuity 



Figure 2. Gyro Noise Model 






Because the four error sources have zero mean value, the expected propagation error is 

zero mean as well. The variance of the error, however, grows with propagation time (t) 
according to the expression 6 w 


°ia( t ) = °Ae(0) + oft + ogt 2 + oj*t 3 /3 


( 1 ) 


The individual error source variances are defined as follows where E is the expectation 
operator and <5 is the Dirac delta function. 


E[b 2 ] = o 2 


( 2 ) 


E[A^(0)] = 4,(0) 


E[n v (t) n v (r)] = o 2 <5(t - r) 


(3) 

(4) 


E[n u (t) n u (r)] = o 2 <5(t - x) (5) 

Standard deviations of the ERBS gyro rates were compiled over the course of the mission 
noise Pr ° Vlde 6Stimate of the 8yro noise - Fi 8 ure 3 shows this history of increasing 

These standard deviations for the sampled angular velocity (a ■ ) give an estimate of the 
float torque noise standard deviation (a v ) that can be used to predict a typical batch 
junction discontinuity. Because the sampled rates are actually averaged over the sampling 
1 1 ”^ rva 'O. t ^ leir var ' ance is smaller than the float torque noise variance by a factor of 


°L, = 05/7 (6) 


Since T is 1 second, the numerical values of a - and <r, are equal, but their units differ 

as thftv snrmlH 


5. EVIDENCE OF GYRO NOISE 

Figure 4 plots daily average junction discontinuities and observation residuals against 
corresponding gyro rate standard deviations. In all cases, the discontinuities and residuals 
increase with increasing noise. 

One thing that requires explanation here is the very steep slope of the yaw and Sun A 
curves. This is due to the low noise on the yaw gyro and the coupling of propagation 
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Figure 3. Gyro Rate Standard Deviations 


errors on the roll and yaw axes. The large roll gyro noise causes yaw as well as roll 
propagation errors. Plotting yaw error against yaw gyro noise alone is therefore mislead- 
ing. A better way to show the dependence might be to combine the roll and yaw into a 
single curve. 

The predicted values of junction discontinuity are based on Equation 1 and use the gyro 
rate standard deviation for <r v • Although inaccurate, the predictions are of the same 
order of magnitude as the observed discontinuities. Also, the rate standard deviation is 
an overestimate of the gyro noise since it includes the actual angular acceleration of the 
spacecraft. With a smaller assumed value for <7 V , the predictions would be closer to the 
observed discontinuities. 
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Figure 4. Discontinuities and Residuals 
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6. THE SIMPLE KALMAN FILTER 


The fact that pitch should be so hard to estimate illustrates one weakness of the batch 
estimator. The reason it is so susceptible to propagation error is that it accepts the 
angular rates almost at face value. After correcting for biases and scale factors, the 
propagation is assumed to be perfect and observations far from the epoch are given equal 
weight with those close to the epoch. 

The Kalman filter differs from the batch estimator in that it estimates the attitude at each 
time point rather than at an epoch time (Reference 7). It also gives greater weight to 
observations made close to the solution time. This seemed to be the answer to the propa- 
gation error/discontinuity problem. Old observations would be forgotten, giving the so u- 
tion more freedom to follow recent observations. Even more, solutions could be made 
continuous at batch boundaries by starting with the previous batch solution and covari- 

ance. 

Although different in approach from the batch estimator, the filter makes many of the 
same computations. For a very simple filter estimating attitude and constant bias, only 
three major modifications are required to the ERBS batch estimator. 

1. Change the partial derivatives 

a. In the batch estimator, the partial derivatives of the current observation 
residual vector Ay^ 


Ay^ = yj (observed) - g (computed) 


( 7 ) 


with respect to the current attitude error vector a, are postmultiplied by the 
epoch-to-current-time propagation matrix 0(tj, to) . This multiplication is 
omitted in the filter 


da 0 


= 0(t j. to) 


dSj 


( 8 ) 


b. Filtered gyro biases bj act only over individual time steps before being 
updated with the bias error vector Oj . Rather than accumulate derivatives 
of the current attitude error with respect to gyro biases from epoch to the 
current time, as in the batch estimator, the filter computes them only for 
the most recent time step. 


d<*j 

dOj 


= tj-i) + 0(t v VO 


daj-i 


tf'Oj. Vi) 


( 9 ) 


104 


The matrix ^/(tj, tj_,) comes from integrating the propagation matrix and 
assumes the bias is constant over this interval. 


V>(tj. <j-i) 



( 10 ) 


c. To simplify notation, the combination of the attitude and gyro bias error 
vectors are referred to as the state Xj 

- [a, ; «j] (11) 


and the partial derivatives of the observation residuals with respect to the 
state is called Fj 


Fj = 



: SC\ m fe") 

' [dxjj 


( 12 ) 


2. Update the state at each time step— The batch estimator weights, transforms, 
and accumulates the observation residuals over the entire batch to give 


FTWY = 


N 


i 


i = 1 


Fj W Ay, 


( 13 ) 


The observation weight matrix W is diagonal with elements equal to the recipro- 
cal of the observation variances 



( 14 ) 


The batch estimator also accumulates the normal matrix, which is the inverse of 
the covariance matrix P. 


i 

i = i 


Fj W Fi + W 0 




( 15 ) 
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The a priori estimate weight matrix W 0 is also diagonal and has as its elements 
the reciprocal of the a priori state variable variances 


w - / 1/o 5> . M (16) 

Wo ■ \ 0 1/<I 

In batch estimation, the epoch state estimate is updated after all the observation 
residuals have been accumulated. 


Ax 0 = (FTWF)" 1 (FTWY) = (P N ) (FTWY) 


(17) 


The filter makes this correction at each time step using the weighted trans- 
formed observation residuals for that instant and the accumulated covariance 
matrix. 


Axj = (Pj) (F/ W Ayj) V*) 

3. Propagate the covariance matrix— The batch estimator propagates the attitude 
and predicts observations at each time step. The filter propagates the co- 
variance matrix as well. This is done by pre- and postmultiplying the 
covariance matrix by the incremental state propagation matrix 0 (tj, tj_i) and its 
transpose. 


Pj = ^(tj. tj-i> Pj-i VO + Q ( 19 ) 

A contribution Q due to the propagation error is also added at this point (Refer- 
ence 3). 


o? (tj - Vi) + \ o 2 (tj - tj .,) 3 - 2 (tj - Vi ) 2 


Q = 


( 20 ) 


-- o 2 (tj - t hl ) 2 


(ft (tj - tj_i) 


The state propagation matrix includes the attitude propagation matrix as its 
upper left block 


<Ktj, tj— i) 


( 


0 (tj. tj-,) | V»( tj. tj-l) 

r - 

0 I I 3 x 3 


) 


( 21 ) 
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The attitude parameterization chosen for the filter is the attitude error vector. This repre- 
sentation fits in with the existing batch differential correction and avoids complications 
from unit length and orthonormality constraints. To update the estimated geocentric 
inertial-to-body attitude matrix A B /i , the error vector is converted to a rotation matrix 
R(§j) which is used to premultiply the old value of A B /i . 


Ab/iOj) = R(aj) Ab/iOj-i) 


( 22 ) 


The gyro bias corrections a } are simply added to the previous bias values 

bj = bj-i + a, (23) 


7. COMPARING THE BATCH ESTIMATOR AND THE FILTER 

After being implemented, the filter was used with low and high noise gyro rate and was 
compared to the results from the batch estimator. For timespans having Sun sensor 
coverage, the attitude is completely observable, and the filter provides smaller root-mean- 
square (rms) observation residuals. This is shown in Table 1. There is no chance of 
junction discontinuities because the filter starts with the last attitude of the preceding 
timespan. 


Table 1. RMS Residuals— Complete Observability 


Observation Residual Variance (deg 2) 


Low Noise 

Roll 

Pitch 

Sun A 

Overall 

Filter 

0.16 

0.30 

0.03 

0.87 

Batch 

0.10 

0.31 

0.01 

0.87 

High Noise 

Roll 

Pitch 

Sun A 

Overall 

Filter 

^ 0.05 

0.13 

mm 

0.17 

Batch 

0.85 

0.73 

0.28 

1.07 
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Over timespans without continuous Sun coverage, the filter can diverge and cause jumps 
even greater than those between batches. An example of such a filter discontinuity is 
provided in Figure 5. There, the yaw diverges until the Sun comes back into view. The 
size of the jump depends strongly on the actual gyro noise and the assumed gyro noise, 
i.e., tuning of the filter. 



Figure 5. Drifting of Filtered Yaw 

As shown in Table 2, the residuals are still much smaller for the filter than for the batch 
estimator. This is in spite of the discontinuity. The reason is that the filter fits all the 
observations better except for the first few following the start of Sun coverage. 
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Table 2. RMS Residuals— Incomplete Observability 


Observation Residual Variance (deg 2) 


Low Noise 

Roll 

Pitch 

Sun A 

Overall 

Filter 

mtm 

0.46 

1.26 

0.40 

Batch 

0.31 

2.58 

0.14 

1.50 

High Noise 

Roll 

Pitch 

Sun A 

Overall 

Filter 

0.04 

0.08 

1.45 

0.19 

Batch 

29.66 

9.67 

0.66 

17.37 


8. DIVERGENCE REMEDIES 


The batch estimator, despite its problems with junction discontinuities and observation 
residuals, does not drift during periods of incomplete attitude observability. This is be- 
cause it has access to data from the entire batch at once. Thus, the batch estimator can 
choose a bias satisfying the observations at both the beginning and the end of the batch. 
The filter, on the other hand, knows only what has been seen up to the current time and 
does not have this advantage. 

To improve filter performance during periods without Sun, the magnetometer can be used 
along with the Sun and Earth sensors. Although it is the least accurate sensor, the mag- 
netometer provides some measure of the yaw at all times. As shown in Figure 6, this 
prevents divergence but does not provide a very accurate definitive attitude solution. 

An alternative to improving the observability is to improve the propagation. To do so, the 
batch estimator can be used to compute a gyro bias over the period lacking Sun coverage. 
This bias can then be used in the filter to reduce the divergence, as shown in Fjlgure 7. 
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Figure 6. Yaw Using Magnetometer 

A variant on this “two-pass” processing is to average the forward-filtered solution with a 
backward-filtered solution obtained by starting from the end of the timespan and letting 
time run backward. This “filter-smoother” can still diverge toward the middle of the Sun 
coverage gap, but does not change abruptly the way a simple filter can. 

The approach that should, perhaps, have been tried first is to tune the gyro noise 
parameter values to reduce the divergence. The values used here for the Q-matrix were 
based on the noise estimated from the raw gyro rates. It was noted that these assumed 
values strongly affect the filter divergence, but no effort was made to find values that 
worked well in all cases. Judging from the satisfactory performance of the ERBS gyro- 
compass, which uses only the Earth sensors and gyro, it should be possible to find such 
values. Perhaps the gyro noise levels used here were too large and so caused the filter to 
forget old observations too quickly. 
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Figure 7. Yaw With Batch Bias 


9. SUMMARY 

This paper uses the ERBS definitive attitude determination system to demonstrate two 
things: first, that gyro noise affects batch attitude accuracy, as seen in the observation 
residuals and batch junction discontinuities; and second, that Kalman filtering can reduce 
those residuals and eliminate the junction discontinuities as well. What was also demon- 
strated, albeit inadvertently, was that during times of incomplete attitude observability, 
the filtered attitude may diverge. Patchwork remedies were tried, such as using additional 
sensors or future observations, but the basic cause of the problem seems to have been the 
“tuning” of the filter. Had a smaller value been assumed for the gyro noise, the attitude 
would not have diverged as quickly as it did. In spite of the present incomplete success, 
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Kalman filtering still promises accuracy equal to or better than that of the batch estimator 
once the tuning process is mastered. 
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